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ABSTRACT 

For the high-frequency peaked BLLac object Mrk421 we study the variation of 
the spectral energy distribution (SED) as a function of source activity, from quiescent 
to active. We use a fully automatized ^-minimization procedure, instead of the "eye- 
ball" procedure more commonly used in the literature, to model nine SED datasets 
with a one-zone Synchrotron-Self-Compton (SSC) model and examine how the model 
parameters vary with source activity. The latter issue can finally be addressed now, 
because simultaneous broad-band SEDs (spanning from optical to VHE photon ener- 
gies) have finally become available. Our results suggest that in Mrk421 the magnetic 
field (B) decreases with source activity, whereas the electron spectrum's break energy 
(7br) and the Doppler factor (5) increase - the other SSC parameters turn out to be un- 
correlated with source activity. In the SSC framework these results are interpreted in 
a picture where the synchrotron power and peak frequency remain constant with vary- 
ing source activity, through a combination of decreasing magnetic field and increasing 
number density of 7 < 7 br electrons: since this leads to an increased electron-photon 
scattering efficiency, the resulting Compton power increases, and so does the total (= 
synchrotron plus Compton) emission. 

Subject headings: BL Lacertae objects: general - BL Lacertae objects: individual (Mrk421) 
- diffuse radiation - gamma rays: galaxies - infrared: diffuse background 



-3- 



1. Introduction 

It's commonly thought that the fueling of supermassive black holes, hosted in the cores 
of most galaxies, by surrounding matter produces the spectacular activity observed in AGNs. 
In some cases 10%) powerful collimated jets shoot out in opposite directions at relativistic 
speeds. The origin of such jets is one of the fundamental open problems in astrophysics. 

If a relativistic jet is viewed at small angle to its axis, the observed emission is amplified 
by relativistic beaming (Doppler boosting and aberration) allowing deep insight into the 
physical conditions and emission processes of relativistic jets. Sources whose boosted jet 
emission dominates the observed emission (blazaro represent a minority among AGN, but 
are the dominant extragalactic source class in 7-rays. Since the jet emission overwhelms all 
other emission from the source, blazars are key sources for studying the physics of relativistic 
extragalactic jets. 

The jets' origin and nature are still unclear. However, it is widely believed that jets are 
low-entropy (kinetic/electromagnetic) flows that dissipate some of their energy in (moving) 
regions associated with internal or external shocks. This highly complex physics is approximated, 
for the purpose of modelling the observed emission, with one or more relativistically moving 



homogeneous pl asma regions (b 



of particles (e.g. iMaraschi et al. 



obs), where radiation is emitted by a non-thermal population 



1992). The high energy emission, with its extremely fast 
and correlated multifrequency variability, indicates that often one single region dominates the 
emission. 

The jet's broad-band (from radio to 7-ray frequencies) spectral energy distribution (SED) 
is a non thermal continuum featuring two broad humps that peak at IR/X-ray and GeV/TeV 



1 Within the blazar class, extreme objects, lacking even signatures of thermal processes usually 
associated with emission lines, are defined as BLLac objects. 
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frequencies and show correlated luminosity and spectral changes. This emission is commonly 
interpreted within a Synchrotron-Self-Compton (SSC) model where the synchrotron and Compton 
peak s are produced by one same time-varying population of particles moving in a magnetic field 



(e.g., 



Tavecchio et al. 



1998, he ncefort T98). The SS C model is perfectly adequate to explain the 



SEDs of BLLac sources (e.g., 



Tavecchio et al. 



20101) 



One important issue that should be addressed, but has not yet been so far because of the lack 
of simultaneous broad-band SEDs, is how the emission changes as a function of the source's 
global level of activity. In particular, given an emission model that fits the data, it is should be 
examined what model parameters are correlated with source activity. In order to investigate 
SEDs at different levels of activity we choose a high-frequency-peaked BLLac (HBL) object, 
i.e. a blazar (i) whose relativistic jet points directly toward the observer, so owing to relativistic 
boosting its SSC emission dominates the source; (ii) whose Compton peak (<; 100 GeV) can be 
detected by Cherenkov telescopes; and (iii) wh ose GeV spectrum can be described as a simple 



power law (unlikely other types of BL Lacs, see 



Abdoetal 



2009J). In addition, such HBL source 



must have several simultaneous SED datasets available. Mrk421 meets these requirements. In 
this paper we study the variation of its SED with source activity, from quiescent to active. 

Another requirement for this kind of study concerns the modeling procedure. We use a 
full-fledged ^-minimization procedure instead of the "eyeball" fits more commonly used in the 
literature. While the latter at most prove the existence of a good solution, by finely exploring the 
parameter space our procedure finds the best solution and also proves such solution to be unique. 

In this paper we investigate the SED of Mrk421 in nine different source states. To this 
aim we fit a one-zone SSC emission model (described in Sect. 2), using a fully automatized 
X 2 -minimization procedure (Sect. 3), to the datasets described in Sect.4. The results are presented 
and discussed in Sect. 5. 
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2. BL Lac SSC emission 



To describe the HBL broad-band emission, we use the one-zone SSC model of T98. This 
has been sho wn to adequately descr ibe broad-band SEDs of most high-frequency-peaked BLLac 



objects (e.g., 



Tavecchio et al 



(ITavecchio et al 



2001 



201 01) and, for a given source, both its ground and excited states 



Tagliaferri et al 



2008|) . The main support for the one-zone model is that 



in most such sources the temporal variability is clearly dominated by one c haracteristic timescale , 



which implies one dominant characteristic size of the emitting region (e.g. jAnderhub et al 



2009). 



Moreover, one of the most convincing evidence favoring the SSC mo del is the strict cor relation 



that usually holds between the X-ray and VHE 7-ray variability (e.g., Fossati et al 
the SSC model the emission in the two bands is produced by the 
and SSC mechanism, respectively), a strict correlation is expecte* 



2008): since in 



the SSC model the emission in the two bands is produced by the same electrons (via synchrotron 

I. 



In our work, for simplicity we used a one-zone SSC model, assuming that the entire SED 
is produced within a single homogeneous region of the jet. As already noted, this class of 
models is generally adequate to reproduce HBL SEDs. However, one-zone models also face 
some problems in explaining some specific features of TeV blazar emission. In particular, 
while very large Doppler factors are often required in on e-zone model, radio VLBI observation s 



hardly de tect superluminal motion at parsec scale (e.g.. (Finer et al 



This led (IGeorganopoulos & Kazanas 



2003; 



Ghisellini et al 



2010; 



Giroletti et al. 



2006)) 



2005) to propose the existence of 



a stru ctured, inhomogeneous and de celerating emitting jet. Inhomogeneous (two-zone) models 



(e.g., (|Ghisellini & Tavecchio 



2008)) have been also invoked to explain the ultra-rap id variability 



occasionally observed in TeV blazars (e.g., (|Aharonian et al. 



2007 



Albert et al. 



2007). 



The emission zone is supposed to be spherical with radius R, in relativistic motion 



2 The rarely occurring "orphan" TeV flares, that are not acco mpanied by variations in the X-ray 
band, may arise from small, low-5, high-density plasma blobs (|Krawczynski et al.ll2004h . 
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with bulk Lorentz factor T at an angle 6 with respect to the line of sight to the observer, so 
that special relativistic effects are cumulatively described by the relativistic Doppler factor, 
5 = [r(l — (3 cos6)Y l . Relativistic electrons with density n c and a tangled magnetic field with 
intensity B homogeneously fill the region. 

The relativistic electrons' spectrum is described by a smoothed broken power-law function of 
the electron Lorentz factor 7, with limits 71 and 72, break at 7t, r and low- and high-energy slopes 
rii and n 2 . This purely phenomenological choice is motivated by the observed shape of the humps 
in the SEDs, well represented by two smoothly joined power laws for the electron distribution. 
Two ingredients are important in shaping the VHE 7-ray part of the spectrum: (i) using the 
Thomson and Klein-Nishina cross sections when, respectively, ^hv < m c c 2 and ^hv > m e c 2 
(with v the frequency of the 'seed' photon), in building the model (T98); and (ii) the correction 



of > 50 GeV data for absorption by the 
photon energy and source dista nce (e.g., 



purpose we use the popular 



xtragalactic Backgroun d Light (EBL), as a function of 



Mankuz hivil et al 



Franceschini et al 



2010 . and references therein): for this 



(2008) EBL model. 



The one-zone SSC model can be fully constrained by using simultaneous multifrequency 
observations (e.g., T98). Of the nine free parameters of the model, six specify the electron energy 
distribution (n e , 71, 7b r , 72, n 1 , n 2 ), and 3 describe the global properties of the emitting region 
(B, R, 5). Ideally, from observations one can derive six observational quantities that are uniquely 
linked to model parameters: the slopes, cti 2 , of the synchrotron bump at photon energies below 
and above the UV/X-ray peak are uniquely connected to rii j2 ; the synchrotron and SSC peak 
frequencies, u SiC , and luminosities, L s C , are linked with B, n c , 5, 7b r ; finally, the minimum 
variability timescale t var provides an estimate of the source size through R^, ct vav 5/ (1 + z). 

To illustrate how important it is to sample the SED around both peaks, let us consider a 
standard situation before Cherenkov telescopes came online, when we would have had only 
knowledge of the UV/X-ray (synchrotron) peak. Clearly this would have given us information on 
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the shape of the electron distribution but would have left all other parameters unconstrained: in 
particular the degeneracy between B and n c - inherent in the synchrotron emissivity - could not 
be lifted without the additional knowledge of the HE/VHE 7-ray (Compton) peak. 

Therefore, only knowledge of observational quantities related to both SED humps enables 
determination of all SSC model parameters. 



3. x 2 minimization 

In this section we discuss the code that we have used to obtain an estimation of the 
characteristic parameters of the SSC model. As we recalled in the previous section the SSC model 
that we are assuming is characterized by nine free parameters, n c , 71, 7 br , 72, nu n 2 , B, R, 5. 
However, in this study we set 71 = 1, which is a widespread assumption in the literature, so 
reducing the number of free parameters to eight. 

The determination of the eight free parameters has been performed by finding their best 
values and uncertainties from a x 2 minimization in which multi-frequency experimental points 
have been fitted to the SSC model d escribed in T98. Minimization has been performed using the 



Levenberg-Marquardt method (see Press et al. 



1994), which is an efficient standard for non-linear 



least- squares minimization that smoothly interpolates between two different minimization 
approaches, namely the inverse Hessian method and the steepest descent method. 

The algorithm starts by making an educated guess for a starting point P in parameter space 
with which the minimization loop is entered; at the same time a (small) constant is defined, 
AxjSjj (where NI stands for 'Negligible Improvement'), which represents an increment in x 2 
small enough that the minimization step can be considered to have achieved no significative 
improvement in moving toward the minimum x 2 '- this will be used as a first criterion for the exit 
condition from the minimization loop. Indeed, to be more confident that the minimization loop 



-8- 



has reached a x 2 value close enough to the absolute minimum, the above condition has to be 
satisfied four times in a row: the number of consecutive times the condition has been satisfied 
is conveniently stored into the integer variable c N j - which is, thus, set to zero at startup. From 
the chosen value of Pq we can compute the associated Xo = X 2 {Po) and enter the minimization 
loop. The Levenberg-Marquardt method will determine the next point in parameter space, P, 
where x 2 Wli l be evaluated. If x 2 {P) > Xo me weight of the steepest-descent method in the 
minimization procedure is increased, the variable cni is set to zero, and we can proceed to the 
next minimization step. If, instead, x 2 {P) < Xo we further check if the decrease in x 2 is smaller 
than or equal to ^Xm- If this is me case, the negligible-improvement counter c NI is increased 
by one: if the resulting value is cni > 4, we think we have a good enough approximation of the 
absolute minimum - and the algorithm ends. If, instead, in the latter test cni < 4, or if in the 
preceding test A% 2 > Axm, men we increase the weight of the inverse Hessian method in the 
minimization procedure, we set x\ equal to the lower value we have just found, and we continue 
the minimization loop. For completeness and illustration , we briefly present the flow-chart of the 
algorithm in Fig. \T\ 

A crucial point in our implementation is that from T98 we can only obtain a numerical 
approximation to the SSC spectrum, in the form of a sampled SED. On the other hand, at each 
step of the loop (see Fig.Q]) the calculation of x 2 requires the evaluation of the SED for all the 
observed frequencies. Usually the model function is known analytically, so these evaluations are a 
straightforward algebraic process. In our case, instead, we know the model function only through 
a numerical sample and it is unlikely that an observed point will be one of the sampled points 
coming from the implementation of the T98 model. Nevertheless it will in general fall between 
two sampled points, which allows us to use interpolation!^ to approximate the value of the SED. 

3 The sampling of the SED function derived from T98 is dense enough, so that, with respect to 
other uncertainties, the one coming from this interpolation is negligible. 
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At the same time, the Levenberg-Marquardt method requires the calculation of the partial 
derivatives of x 2 with respect to the eight fitted SSC parameters. Contrary to the usual case, in 
which from the knowledge of the model function all derivatives can be obtained analytically, 
in our case they have also been obtained numerically by evaluating the incremental ratio of the 
X 2 with respect to a sufficiently small, dynamically adjusted increment of each parameter. This 
method could have introduced a potential inefficiency in the computation, due to the recurrent 
need to evaluate the SED at many, slightly different points in parameter space, this being the most 
demanding operation in terms of CPU time. For this reason we set up an algorithm to minimize 
the number of calls to T98 across different iterations. 



4. Datasets 



From the literature we then select nine SED datasets corresponding to different emission 
states (low to high) of the HBL source Mrk421. 



State 1 and state 2 (|Acciari et al. 



2009) multi-wavelength campaigns were triggered by 



a major outburst in April 2006 that was detected by the Whipple 10 m telescope. A prompt 
campaign was not possible because of visibility constraints on XMM-Newton. So simultaneous 
multi-wavelength observations took place during the decaying phase of the burst. The optical/UV 
and X-ray observations were carried out using XMM-Newton's optical monitor (OM) and 
EPIC-pn detector, respectively. The MAGIC and Whipple telescopes were used for the VHE 7-ray 
observations. State 1 strictly simultaneous observations lasted ~4hrs for state 1, and more than 
3 hrs for state 2. 



State 3 (Rebillot et al. 



2006) reports the multi-wavelength observations during December 



2002 and January 2003. The campaign was initiated by X-tay and VHE flares detected by the 
All Sky Monitor (ASM) of the Rossi X-ray Timing Explorer (RXTE) and the 10 m Whipple 
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telescope. Whipple and HEGRA-CT1 were used for the VHE observations during the campaign. 
Even though Whipple observed the source from Dec. 4, 2002 to Jan. 15, 2003, and HEGRA-CT1 
on Nov.3-Dec.12, 2002, only the data taken during nights with simultaneous X-ray observations 
are used in this paper to construct the SED. Optical flux is the average flux of the data obtained 
from Boltwood observatory optical telescope, KVA telescope, and WIYN telescope during the 
campaign period. 



State 4 and state 9 (|Blazejowski et al 



2005|) observations were taken from a comparatively 



longer time campaign in 2003 and 2004. The X-ray flux obtained from 7?XTE were grouped 
into low-, medium- and high-flux groups. For each X-ray observation in a given group, Whipple 
VHE 7-ray data that had been observed within an hour of the X-ray data was selected. State 4 (i.e., 
medium-flux) was observed between March 8 and May 3, 2003; whereas state 9 (i.e., high-flux) 
was observed on April 16-20, 2004. Optical datasets obtained with Whipple Observatory's 1.2 m 
telescope and Boltwood Observatory's 0.4 m telescope were also selected based on the same 
grouping method. The optical data measured during the whole campaign were not simultaneous 
with the other multi- wavelength data: however, the optical flux was not found to vary significantly 
during the campaign, so its highest and lowest values are taken to be reliable proxies of the actual 
values. 



State 5 and state 7 (|Fossati et al.1 120081) data were taken on March 18-25, 2001 during a 
multi-wavelength campaign. State 7 denotes the peak of the March 19 flare, whereas state 5 
denotes a post-flare state on March 22 and 23. In both cases the X-ray and VHE 7-ray data were 
obtained with, respectively, 7?XTE and the Whipple telescope. The lowest and highest optical 
fluxes obtained during the whole campaign with the 1.2 m Harvard-Smithsonian telescope on 
Mt. Hopkins were used in the SEDs for states 5 and 7, respectively. 



State 6 ( Acciari et al 



2009) observations were taken, using the same optical and X-ray 



instruments as states 1 and 2, during a decaying phase of an outburst in May 2008. The VHE 7-ray 
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data were taken with VERITAS. There are ~2.5 hours of strictly simultaneous data. 



State 8 (Donnarumma et al 



2009) data were taken during a multi-wavelength campaign on 



June 6, 2008: VERITAS, RXTE and Swift/BAT, and WEBT provided the VHE 7 -ray, X-ray, and 
optical data, respectively. 



5. Results and discussion 



To each of the datasets we apply our ^-minimization procedure (see Fig. \\}. The best-fit 
SSC models are plotted alongside the SED data in Fig. |2] 

Obtaining truly best-fit SSC models of simultaneous blazar SEDs is crucial to measure 
the SSC parameters describing the emitting region. The obtained optimal model is proven to 
be unique because the SSC manifold is thoroughly searched for the absolute x 2 minimum. 
Furthermore, the very nature of our procedure ensures that there is no obvious bias affecting the 
resulting best-fit SSC parameters. 



20_10j), there may be caveats 



Nevertheless, as it has also been discussed (lAndrae et al. 
related with the xl fitting, especially when applied to non-linear models such as the present one. 
For this reason, it is important to try to understand th e goodness of the fi t with methods other 



than the value of xl (reported in tabled). Following (|Andrae et al 



2010|) . we have applied the 



Kolmogorov-Smimov (KS) test for normality of the residuals of all our SED fits. A standard 
application of the test shows that in all cases the residuals are not normally distributed: the KS 
test thus fails at the 5% significance level. It is, of course, crucial to understand the reason for 
this behavior of our SSC fits with respect to the KS test. Let us start with some general remarks 
about the modelling of blazar emission and their observations. First, the one zone SSC model 
contains two distinct physical processes in one same region, i.e. synchrotron emission and its 
Compton up-scattered counterpart, that manifest themselves as essentially separate components 
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at very different energies; on the other hand, additional subtle effects may enter the modelling of 
blazar emission, so that the SSC model may only be an approximation to the real thing. (A more 
refined KS analysis suggests exactly this, see below). Second, our blazar datasets do cover the (far 
apart) energy ranges spanned by, respectively, the synchrotron and Compton emission: but they 
markedly differ in these two spectral regions, in that the uncertainties associated with VHE data 
are much larger than those associated with the optical and X-ray data. 

Both these observations suggest us and give us the possibility of a slightly different approach 
to the problem of the statistical significance of the fits, which will turn out to be quite enlightening. 
We will refer to this other approach as the piecewise KS test: it consists in applying the KS test, 
separately to low energy and high energy data (see Appendix). The main motivation behind this 
idea is, as discussed above, the marked difference, from both the physical and the data-quality 
point of view, of the lower and higher energy ranges. Surprisingly, if for each SED we separately 
check the low- and high-energy residuals for normality, the KS test always confirms their 
(separate) normality at the 5% confidence level. On the technical side we took all the necessary 
precautions to improve reliability: the null hypothesis statistics have been obtained from a 
Montecarlo simulation of 100,000 datasets having the same dimension as the residuals datasets. 
Critical values to test normality of the residuals have been obtained from these numerically 
determined statistics and, normality holds in all that cases we have considered. Hence, it is 
unlikely to be a coincidence, and this calls for an explanation. Clearly, the fact that the piecewise 
KS test is not able to reject the normality of the low-/high-energy residuals separately, means that 
the fitted SEDs can be considered as a reasonable models separately at low and high energies. Of 
course, because of the quality of especially high energy data, the uncertainties on the parameters 
of the fit are sometimes quite large, a fact that unfortunately can not be avoided, despite the fact 
that the quality of the datasets that we are using is certainly above average among those that are 
available. At the same time, the failure for normality of the residuals on the standard KS test, in 
which low- and high-energy residuals are considered simultaneously, suggests that existing data 
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may require our adopted SSC model to be improved - perhaps by taking into account higher-order 
effects (e.g., a multiple-break or curved electron spectrum). So, and with more and better VHE 
data available, it could be possible to reduce the uncertainties in the parameters and to obtain a 
higher overall significance in the fit. 

Given the above, the results found with our fits are the best possible in the framework of 
existing datasets and models of blazar emission: at least, although preliminary, they have a 
quantitative and clear statistical meaning. Therefore, we think it is useful to examine some of their 
possible consequences. We'll do so in what follows. 

In Tables [B HI we report the best-fit SSC parameters. Source activity (measured as the total 
luminosity of the best-fit SSC model) appears to be correlated with B, 7 br , and S (see Fig.[3]-to/?) 
- and to be uncorrected with the remaining SSC parameters. The bolometric luminosity used 
in these plots has been obtained directly from the fitted SED. In more detail, after determining 
the numerical approximation to the SED \og[uF(u)}, the parameters being fixed at their best 
values obtained with the previously described minimization procedure, we have performed 
L = T' /max vF(v)dv, with z/ min , z/ max set at 2.5 decades, respectively, below the synchrotron peak 
and above the Compton peak. In this way we make sure to perform the integral over all the 
relevant frequencies in a way that is independent from any in location of these peaks. 

We then searched the data plotted in the top row of Fig. |2]for possible correlations. The 
linear-correlation coefficients turn out to be 0.67, 0.64 and 0.54, which confirm linear correlations 
with confidence levels of 4.8%, 6.3% and 13.3%, respectively. As an additional test (given the 



relatively 
residuals^. 



n 



ow statistics of our datasets), we checked that the KS test confirms normality of the fit 



4 The same approach involving a Montecarlo generated empirical distribution for the null- 
hypothesis described just above has been used also in all these cases. 
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All the parameters derived through our automatic fitting procedu re are within the rang e 



of SSC parameters found in the literature for HB 



2010 



Tagliaferri et al 



2008 



(Bednarek & Protheroe 



Konopelko et al 



2003; 



Celotti & Ghise 



1997 



Tavecchio et al 



Fossati et al 



mi 



s in general (e.g. (ITavecchio et al. 



2001 



1998 



2008) ) and for Mrk421 in particular (e 



Maraschi et al. 



1999; 



Ghisellini et al. 



2002; 



20081 )). In particular, large Doppler fact ors such as those 



derived in our extreme c ases, 5 > 50, have been oc casionally derived (e.g ., (|Konopelko et al. 



2003 



Fossati et al. 



2008); see also the discussion in ( Ghisellini et al 



2005)) 



As is seen from Fig.|3]-top, 7b r and B are correlated, respectively, directly (left) and inversely 
(right) with L. This may be explained as follows. An increase of 7 br implies an effective increase 
of the energy of most electrons (or, equivalently, of the density of 7 < 7b r electrons). To keep the 
synchrotron power and peak roughly constant (within a factor of 3; see Fig. |2]>, B must decrease. 
This improves the photon-electron scattering efficiency, and the Compton power increases. The 
total (i.e., synchrotron plus Compton) luminosity will be higher. So a higher 7 br implies a lower 
B and a higher emission state. The 5-L correlation (Fig. ^-top-middle) results from combining 
the B—6 inverse correlation (Fig. ^-bottom-right; see below) and the B-L anticorrelation. 

A deeper insight on emission physics can be reached plotting the three L-dependent 
parameters one versus the other (Fig. ^bottom). The B-q hr anticorrelation, with AlogB ~ 
—2 Alog7br (Fig. ^-bottom-left), derives from the synchrotron peak, u s oc B7 2 , staying roughly 
constant (see Fig. |2]). For fixed synchrotron and Compton peak frequencies in a relativistically 
beamed emission, the B-5 relati on is predicted to be i nverse in the Thompson limit and direct 



in the Klein-Nishina limit (e.g., 



Tavecchio et al 



1998): because u s , v c do not greatly vary 



from state to state in our data (see Fig. |2]), the correlation in Fig. ^-bottom-right suggests that 
the Compton emission of Mrk421 is always in the Thompson limit. The <5-7b r correlation 
(Fig. 3-bottom-middle) results as a corollary of the condition of constant v B , u c emitted by a 
plasma in bulk relativistic motion toward the observer. 
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Our fits show clear trends among some of the basic physical quantities of the emitting region, 
the magnetic field, the electron Lorentz factor at the spectral break, and the Doppler factor (see 
Fig. [3]). In particular, B and 7b re ak follow a relation B oc %^ eak , wn il e B and 5 are approximately 
related by B oc 5~ 2 . 

Rather interestingly, the relation connecting B a nd 7 is naturally expected within the context 



of the simplest electron acceleration scenarios (e.g., (|Henri et al. 



1999)). In this framework, the 



typical acceleration timescale, t acc , is proportional to the gyroradius: t acc (7) oc r L /c, where 
tl = , ym e c/ (eB) is the Larmor radius. On the other hand, acceleration competes with radiative 
(synchrotron and IC) cooling. In Mrk421, characterized by comparable power in the synchrotron 
and IC components, we can assume that £0001(7) ~ ^ s Y n °^ 1/7-S- The maximum energy reached 
by the electrons is determined by setting these two timescales equal, i.e. 

£acc(7max) £cool(7max) ^ 7max OC B ^ } (1) 

in agreement with the relation derived by our fit. It is therefore tempting to associate 7break to 7 max 
and explain the B oc %^ eak relation as resulting from the acceleration/cooling competition. 

If the above inference is correct, one can explain the variations of 7b rea k as simply reflecting 
the variations of B in the acceleration region. In turn, the variations of B could be associated 
either to changes in the global quantities related to the jet flow or to some local process in the jet 
(i.e. dissipation of magnetic energy through reconnection). The second relation mentioned above, 
i.e. 5 vs. B, seems to point to the former possibility. First of all let us assume 5 ~ T (since we are 
probably observing the Mrk421 jet at a small angle w.r.t. the line of sight). A general result of jet 
acceleration models is that, during the acceleration phase, the jet has a parabolic shape, R oc d 1 ! 2 



(with d the distanc e from the central b 



as T oc d 1 ' 2 (e.g., dviahakis & Konigl 



ack hole), and the bulk 



2004; 



Komissarov et al 



.orentz factor T increases with d 



20071 ). On the other hand, if the 



magnetic flux is conserved, B oc R 2 oc d 1 and thus we found B oc T 2 . Albeit somewhat 
speculative, these arguments suggest that the trends displayed in Fig. [3] are naturally expected in 
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the general framework of jet acceleration. 

The correlations in Fig. ^bottom seem to be tighter than those in Fig. |3]-top. The larger 
scatter affecting the latter owes to the fact that the electron density, n e , that also enters the 
definition of SSC luminosity, shows no correlation with B, 7 br , and 8 - hence it slightly blurs the 
latter's plots with luminosity. 

One further note concerns error bars. Our code returns 1-er error bars. To our best knowledge, 
this is the first time that formal errors of SED fits are obtained in a rigorous way. As an example 
of the soundness of the method, note that the obtained values of 5 are affected by the largest errors 
when the distribution of the VHE data points is most irregular (e.g. states 4, 7). 

We notice that the variability of Mrk421 markedly differs from that of (e.g.) the other nearby 
HBL source, Mrk501. The extremely bursting state of 1997 showed a shift of v s and v c by two 
and o ne orders of mag nitude, respectively, suggesting a Klein-Nishina regime for the Compton 
peak (|Pian et al.lll998|) . Based on the data analyzed in this paper, Mrk421 displays (within 
our observational memory) a completely different variability pattern. However, one important 
similarity may hold between the two sources: based on eyeball-fit analysi s of Mrk 50 l's SEP in 
different emission states using a SSC model similar to the one used here, lAcciari et al.l (|2010l) 
suggest that 7 br does vary with L. If this correlation is generally true in blazars, the implication 
is that particle acceleration, providing fresh high-energy electrons within the blob, must be one 
defining characteristic of excited source states. 
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anonymous referee for useful comments and suggestions. One of us (SA) acknowledges partial 
support from the long-term Workshop on Gravity and Cosmology (GC2010: YITP-T- 10-01) at 
the Yukawa Institute, Kyoto University, during the early stages of this work, and warmly thanks 
P. Creminelli and S. Sonego for insightful discussions on some topics touched upon in this paper. 
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START 

assign Po, initial point 
in SSC parameter space, 
and Axjjjj; set cni = 



cni = 

xl = x 2 {P) 



compute X 2 (Po) 
X 2 o = X 2 (Po) 



determine P, next point 
in SSC parameter space 



output values of x 2 , 
SSC parameters and 
associated uncertainties 
END 



± 



Xo = 


x 2 (P) 


increase 


weight of 


inverse 


Hessian 



NO 



NO 



YES 

cni < 4 



increase weight of 
steepest descent 



x 2 (P) > x 2 



; compute \ ( P) j- 



x 2 (P) < xl 



X 2 (P)-Xo<^Xm : 
YES 



J 



cni = cni + 1 



Fig. 1 . — Flow chart of the minimization code. 



Source 


B 


R 


5 






[gauss] 


[cm] 






State 1 


(9 ±3) x 1CT 1 


(9 ± 4) x 10 14 


(2 ±0.5) x 10 1 


0.84 


State 2 


(8 ±6) x 1CT 1 


(8 ± 4) x 10 14 


(2.7 ± 1.1) x 10 1 


1.86 


State 3 


(6 ±6) x 1(T 2 


(2.0 ±1.5) x 10 15 


(1.0 ±0.5) x 10 2 


0.91 


State 4 


(1.21 ±0.16) x 1CT 1 


(1.1 ± 1.3) x 10 15 


(8 ±6) x 10 1 


0.89 


State 5 


(1.9 ±1.3) x 1CT 1 


(10 ±4) x 10 14 


(7 ±5) x 10 1 


0.67 


State 6 


1.0 ±0.7 


(6 ± 3) x 10 14 


(2.8 ± 1.1) x 10 1 


1.39 


State 7 


(4 ±3) x 10~ 2 


(2 ± 5) x 10 15 


(8 ±7) x 10 1 


1.61 


State 8 


(6 ±3) x 10~ 2 


(2 ±1.8) x 10 15 


(1.1 ±0.4) x 10 2 


0.60 


State 9 


(4 ±3) x 10" 2 


(2 ±4) x 10 15 


(1.2 ± 1.0) x 10 2 


0.85 



Table 1: Best-fit single-zone SSC model parameters for the nine datasets of Mrk421. States are 
named as in Fig. [21 
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Fig. 2. — Best-fit one-zone SSC models for nine data sets referring to different emission lev- 



els of the HBL source Mrk421. Source states are ordered 



data have been obtained as follows: state 1 (Acciari et aL 



sta te 3 IR ebillot et al. 



6 (Acciari et al 



20061), state 4 (Blazeiow ski etal. 



2009), state 7 (Fossati et al 



(IBlazejowski et al.ll2005|) 



jy in creasing model luminosity an d 



2009), state 2 (Acciari et al. 



2005h . state 5 (Fossati et al. 



2008), state 8 (|Donnarumma et al 



2009), 



2008), state 



2009), state 9 
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4.5 5 
L °g ?break 



L °g 7 break 




Fig. 3. — Top. Variations of the SSC parameters magnetic field B, Lorentz factor 5, and 7 br as a 
function of the model's bolometric luminosity. The other SSC parameters that are left to vary (R, 
n e , 72, Tii, fi2 only show a scatter plot with the luminosity. Bottom. Correlations between B, 5, and 



7br 
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APPENDIX 

A TOY MODEL MIMICKING THE PIECEWISE KS TEST RESULTS FOR THE SED'S 

In this appendix we reproduce the outcome of the piecewise KS approach to goodness of fit 
using a simulated toy model. Let us consider the function p(x; a, b, c) := ax 4 + bx 2 + a By using 
a pseudo-random generator, that produces uniformly distributed pseudo-random numbers in the 
interval [0, 1] - which we will henceforth call rnd -, we generated (x i: yi)i=i,...,5 pairs of points and 
(u>i)i=i,...,5 numbers, where Xi = — 2.3+0. l-(i+0.1-rnd), yi = p(— 2.3+0. —1,4, 0)+0.2-rnd 
and wi = a^ 2 with Oi = 0.01 ■ (rnd + 1). These can be considered as 5 measurements around the 
x = — 2 local maximum of the function — x 4 + 4a; 2 with small uncertainties ov We then generated 
another 5 pairs of points (x i: yj)i=6,...,io an d {wi)i=&,...,io numbers with = 1.7+0.1-(i+0.1-rnd), 
Di = p(1.8 + 0.1 • i; -2, 16, -14) + 0.3 • rnd, w t = a^ 2 and, finally, a t = 0.1 • (rnd + 1), which 
again can be considered as 5 slightly randomized measurement around a maximum, but now the 
maximum of the function — 2x 4 + 8x 2 + 1 and with much larger uncertainties that the first five 
ones. 

By minimizing x 2 we can fit the above set of 10 randomly generated data points versus the 
function p(x; a, b, c). After obtaining the desired fitted parameters we calculate the residuals and, 
using the KS test, check if they are normally distributed. The KS test rejects normal distribution 
of the entire set of residuals at the 5% significance level. On the other hand, if we separately apply 
the KS test to the first, i.e. {i — 1, . . . , 5}, and the second, {i — 6, . . . , 10}, subset of residuals, 
the test confirms that the residuals are normally distributed. This exactly replicates the behavior 
we obtained in the SED fits. 
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Source 


71 c 


7br 


7max 


ni 


n 2 




[cm -3 ] 










State 1 


(1.3 ±1.5) x 10 3 


(2.6 ±0.9) x 10 4 


(1.05 ±0.18) x 10 7 


1.49 ±0.19 


3.77 ±0.11 


State 2 


(1 ±3) x 10 3 


(2.4 ±0.9) x 10 4 


(4.1 ± 1.1) x 10 6 


1.5 ±0.3 


3.62 ±0.14 


State 3 


(5 ± 5) x 10 3 


(7 ±3) x 10 4 


(7 ±5) x 10 7 


2.05 ±0.10 


4.8 ±0.3 


State 4 


(2 ± 5) x 10 3 


(4 ± 2) x 10 4 


(8.2 ±1.7) x 10 6 


1.8 ±0.3 


4.11 ±0.13 


State 5 


(2 ± 5) x 10 3 


(4.5 ±1.9) x 10 4 


(2.4 ±0.3) x 10 7 


1.7±0.3 


4.3 ±0.180 


State 6 


(4 ± 4) x 10 3 


(1.9 ±0.6) x 10 4 


(1.8 ±0.4) x 10 6 


1.54 ±0.11 


4.37 ±0.09 


State 7 


(1 ± 7) x 10 3 


(8 ± 6) x 10 4 


(7 ±2) x 10 6 


1.7 ±0.4 


4.23 ±0.20 


State 8 


(4 ±9) x 10 1 


(5 ±2) x 10 4 


(1.6 ±0.4) x 10 7 


1.5 ±0.2 


4.22 ±0.14 


State 9 


(1 ± 7) x 10 2 


(8 ± 9) x 10 4 


(1.1 ±0.4) x 10 7 


1.6 ±0.5 


3.9 ±0.2 



Table 2: Best-fit single-zone SSC model parameters for the nine datasets of Mrk421. Numbering 
convention as in Fig. [2] 
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